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ABSTRACT 

» / 

# ^ ' 

Thitpaper presents the efforts of Mechanical Dynamics, Inc. in obtaining a general formulation for flexible 
bodies in a multibody setting, The efforts being supported by MDI. both in house and externally are 
summarized. The feasibility of using lumped mass approaches to modeling flexibility in a multibody dynamics 
context is examined. The kinematics and kinetics for a simple system consisting of two rigid bodies connected 
together by an elastic beam are developed in detail. Accuracy, efficiency and ease of use using this approach are 
some of the issues that are then looked at. 

The formulation is then generalized to a "superelement" containing several nodes and connecting several bodies. 
Superelement kinematics and kinetics equations are developed. 

The feasibility and effectiveness of the method is illustrated by the use of some examples illustrating 
phenomena common in the context of spacecraft motions. 
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PROFILE OF MECHANICAL DYNAMICS, INC. 

COMPANY BACKGROUND 

• HISTORY 

• PRODUCTS & SERVICES 

• CUSTOMERS 

CURRENT PRODUCTS 

• ADAMS 

• ADAMS /MODAL 

• POST PROCESSORS 

SERVICES 

• CONSULTING 

• TRAINING 

• HOTLINE 

AVAILIBILITY OF PRODUCTS 
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MDI EFFORTS IN FLEXIBILITY 

INTERNAL R&D 

• LUMPED MASS APPROACHES TO FLEXIBILITY 

EXTERNAL R&D 

• UNIVERSITY FUNDED RESEARCH IN MODAL APPROACHES 

INTERFACE TO FEA PROGRAMS 

• NASTRAN 

• ANSYS 
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LUMPED MASS APPROACH TO FLEXIBILITY 
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BEAM DEFORMATION DURING MOTION 
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TRANSLATIONAL DISPLACEMENT COMPUTATION 



G 


S 12 - 1*2 + r 2 ' ^1 • Tj 

i£n = A 10 [ B 2 +A 2G t2 • Bi ■ A 1G c 1 ] 
I A 2*2 = 1 £ 12 ■ 1 L 
lL = {L # 0 }T 
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ANGULAR DISPLACEMENT COMPUTATION 
SPACE 1-2-3 ANGLES ARE USED FOR MEASURING ANGLES 



c 2 c 3 

Si S 2 C 3 - S 3 C x 

Ci s 2 C 3 + S 3 Si 

> 

N 

II 

c 2 s 3 

Si S 2 S 3 - C 3 Si 

Ci s 2 s 3 + C 3 Si 


■ s 2 

S1C2 

Cl c 2 

B 2 = SIN - a 31 ) 



IF fi 2 * 

n/2 THEN 



= ARCTAN 2 ( a 32 , a 33 ) 

B3 — ARCTAN2 ( a 12 + a 13 f a 13 - a 22 ) ■ ®i 
ELSE IF B 2 = n /2 THEN 

fi 3 = ARCTAN2 ( 2 21 ( ^ ji ) 
fi 3 = ARCTAN2(a 12> a 13 ) + B 3 

li 2'2 = t ®1 ^2 ^3 
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VELOCITY COMPUTATION 



m 2 , 2 = «3)v 21 . Wj x s 21 

il 2 ' 2 = A u [ A 1G (Ki - ) ■ A JI ffl-in 

- ( A JI ti + A jg Ki - A jg Ri) ] 

CO12 = - CO j 

lfiJl2= A u [ A^SSi - 22 j ] 
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FORCE COMPUTATION 

FORCES AT ORIGIN OF COORDINATE SYSTEM ON REF, FRAME 1 

1 E 2 = - [ Kll lA2'2 + &12 ^ 12 -J 

- [ Qu Il2'2 + £l2 2^121 

1 1 2 = ’ [ K 21 1 A 2’2 + K 22 li 12 ] 

- [ £ 21 1 1 2'2 + £ 22 1 ^ 12 1 

K IS THE STANDARD MATRIX FOUND IN ANY STRUCTURAL ANALYSIS TEXT 


FORCES AT ORIGIN OF COORDINATE SYSTEM ON REF. FRAME 2 
SINCE THE BEAM IS MASSLESS, APPLYING LAWS OF EQUILIBRIUM : 

Fj + F 2 =0 

1 E 1 = - 1 E 2 

Xj + X2 + S Y2 X F2 = 0 

lll = ‘[ll2 + 1 S~12 1 E 2 1 
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ACCURACY OF METHOD 


• DIRECTLY RELATED TO DEGREE OF DISCRETIZATION 


• METHOD DOES NOT YIELD WRONG ANSWERS 


• DEGREE OF DISCRETIZATION DEPENDENT ON FREQUENCY 

CONTENT DESIRED. ADAMS/MODAL WELL COMPUTE 
EIGENVALUES AND EIGENVECTORS FOR ANY ADAMS MODEL. 
CAN ANIMATE LINEAR MODEL USING SELECTED SET OF MODE 
SHAPES AND FREQUENCIES. 
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FLEXIBLE BODY DYNAMICS 

Benchmark Problem #1: Uniform Cantilever Beam Spin-up 


The inset of Fig. 1 shows a uniform, homogeneous, cantilever beam supported by a circular hub of 
radius r. At time t = 0, the system is at rest in a Newtonian reference frame and the beam is undeforme ^ 
Subsequent to this initial time, the hub is made to rotate about a vertical axis X - X, passing through th 
center of the hub, in such a way that ft, the angular speed of the hub, is given by 


0(0 


_ f (2/5) [< - (7.5/ x) sin(*t/7.5)] rad/sec 0 < < < 15 sec 


\6 rad/i 


sec 


t > 15 sec 


which represents a smooth transition from zero hub motion to a constant angular speed of 6 rad/sec. The 
beam has a length L, Young’s modulus E, shear modulus G, mass per unit length p, and a circular cross- 
section of area A and area moment of inertia I. 

The solid line in the figure below shows the time history of the displacement of the beam tip. in the 
plane of rotation, relative to a line fixed in the hub and originally parallel to the centroidal axis of the bearm 
This result was obtained using the theory and algorithm presented m Refs.[l] and [2] with three assume 
modes and the following parameter values 

r = 0m p= 1.2 kg/m 
L = 10 m A - 4 x 10“ 4 rn 2 

£ = 7xl0 10 N/m 2 G = 3 x lO 10 N/m 2 7 = 2xl0" 7 m 4 

All external forces were neglected and the assumed modal functions were chosen to be equal to the first three 
eigenfunctions of an identical uniform cantilever beam with Us root fixed in a Newtonian reference rame. 
The numerical integration was carried out using a 4‘* - 5‘ h order, variable step-size Runge-Kutta-Merson 
method with a print step and initial time step of .03 seconds and an error tolerance of l x 10 . The dashed 

line result was produced with an algorithm based on the assumed-mode formulation utilized in most flexible 
multibody programs. This result was verified by Fidelis Eke [(818) 354-2916] at Jet Propuls, on Labs usmg 

DISCOS 



[1] Ryan, R.R., “Flexibility Modeling Methods in Multibody Dynamics,” Ph D. dissertation, Stanford 
University, Available from University Microfilms, Ann Arbor, Michigan, 1986. 

[2] Kane, T.R., Ryan. R.R., and Banerjee, A.K., “Dynamics of a Cantilever Beam Attached to a Moving 
Base.” Journal of Guidance. Control, and Dynamics , Vol.10, No.2. March-Apnl, 1987. 


2 55 








UNIFORM BEA 
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EFFICIENCY OF METHOD 


• THIS METHOD IS USABLE FOR SMALL TO MEDIUM SIZE PROBLEMS 
( MEDIUM = 300 RIGID AND FLEXIBLE DOF ) 

• FOR LARGER PROBLEMS IT MAY PROVE TO BE MORE CPU INTENSIVE 
THAN DESIRABLE. 


• THE CPU TIME TAKEN FOR A SIMULATION IS LINEARLY 

PROPORTIONAL TO THE NUMBER OF FLEXIBLE BEAMS IN 
THE SYSTEM 
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Number of Beams Us. CPU Time 
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Number of Beams 



Number of Equations Us. CPU 
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EASE OF USE 


• THE RESULTING PROGRAM IS EXTREMELY EASY TO USE. 

• USERS DO NOT NEED STRONG FEA BACKGROUND TO CREATE 
MODELS OF STRUCTURES 

• RECOGNITION AND SELECTION OF PROPER MODES IN AN ART. 
THE RESULTS ARE ONLY AS GOOD AS THE SELECTED MODES. 
DIFFICULTY ALLEVIATED IN THIS APPROACH. 
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GENERALIZATION TO SUPERELEMENTS 



NODES 

• LOCATION OF NODES 2, 3, 4 WRT. TO A KNOWN 
REFERENCE FRAME 

• MASS AND INERTIA PROPERTIES FOR EACH NODE 
OBTAINED FROM MASS MATRIX. 

LOCATION AND ORIENTATION OF COORDINATE SYSTEMS 

2 ON PART 2 2’ ON PART 1 

3 ON PART 3 3' ON PART 1 

4 ON PART 4 4' ON PART 1 
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ASSEMBLY 

♦ PRE-TENSION AND INITIAL DISPLACEMENTS 
AT CONNECTION POINTS 2-2', 3-3', 4-4' 

FLEXIBILITY PROPERTIES 

• STIFFNESS MATRIX 
. DAMPING MATRIX 

DISPLACEMENT COORDINATES 

] T 

__ _ _ _ _ _ _ . 

H a [ H 2'2 3'3 4'4 1 T 

VELOCITY COORDINATES 

1 X = [ i X 2'2 1 X 3'3 1 ¥ 4'4 1 T 

1 = [ i 2'2 1 3'3 1 ® 4'4 J 

FORCE DEFINITION AT COORDINATE SYSTEMS 2, 3, 

lE = [ ] E 2 1E3 1E4 ] T 
iX = [1X2 1I3 lid 7 
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FORCE COMPUTATION AT COORDINATE SYSTEMS 2, 3, 4 

lE = - 1 K 11 1 A + & 12 fi] 

■[ £n iY +£i2 ifiJ 1 +Eo 

ll = ■ [ K21 i A +K22 fi ] 

"I El 21 i¥- + £ 22 i& 1 +Io 

FORCE AT COORDINATE SYSTEM 1 
F| + F2+ F3+ F4 = 0 

1E1 = • t iE 2 + 1E3 + 1E4] 

Ti + T 2 + T3 + T 4 + S 12 x F 2 +Sj 3 x Fjf s 14 x F 
1X1 = - [ 1E2 + 1E3 + ie 4 ] 

- [lS'12 1E2 + I5‘l3 1E3 + lS'14 1E4] 
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